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A Wigner-Poisson kinetic equation describing charge transport in doped semiconductor 
superlattices is proposed. Electrons are supposed to occupy the lowest miniband, ex- 
change of lateral momentum is ignored and the electron-electron interaction is treated 
in the Hartree approximation. There are elastic collisions with impurities and inelas- 
tic collisions with phonons, imperfections, etc. The latter are described by a modified 
BGK (Bhatnagar-Gross-Krook) collision model that allows for energy dissipation while 
yielding charge continuity. In the hyperbolic limit, nonlocal drift-diffusion equations 
are derived systematically from the kinetic Wigner-Poisson-BGK system by means of 
the Chapman-Enskog method. The nonlocality of the original quantum kinetic model 
equations implies that the derived drift-diffusion equations contain spatial averages over 
one or more superlattice periods. Numerical solutions of the latter equations show self- 
sustained oscillations of the current through a voltage biased superlattice, in agreement 
with known experiments. 

1. Introduction 

Nonlinear charge transport in semiconductor superlattices has been widely stud- 
ied in the last decade ' 2 . A superlattice (SL) is a convenient approximation to 
a quasi-one-dimensional crystal that was originally proposed by Esaki and Tsu to 
observe Bloch oscillations, i.e., the periodic coherent motion of electrons in a mini- 
band in the presence of an applied electric field. When the materials were grown, 
many interesting nonlinear phenomena were observed, including self-oscillations of 
the current through the SL due to motion of electric field pulses, multistability of 
stationary solutions, and so on 1,2 . In the effort to describe quantitatively trans- 
port in these materials, a large number of quantum and semiclassical equations, 
hydrodynamic and drift-diffusion models have been proposed. However, systematic 
derivations of reduced balance equations (hydrodynamic or drift-diffusion equa- 
tions) are scarce. Hydrodynamic equations have been derived by Lei and coworkers 
3,4 and their numerical solutions describe self-oscillations of the current in different 
regimes 5 . They write operator Heisenberg equations, and average them to obtain 
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a hierarchy of equations for moments such as electron density, current and energy 
These equations are then closed by using a Fermi-Dirac local equilibrium ansatz, 
which clearly does not amount to a systematic procedure, nor gives a clue to the 
validity of the resulting hydrodynamic equations. A similar closure procedure has 
been recently developed for the case of quantum particles in an arbitrary external 
three-dimensional potential by Degond and Ringhofer 6 , who have used the maxi- 
mum entropy principle to select their local equilibrium. These authors employ an 
operator form of a local Maxwellian in their examples, which therefore do not obey 
the Pauli exclusion principle. Earlier, Huang and Wu 7 found a local Fermi-Dirac 
distribution function as a basis of their hydrodynamic equations by using a different 
definition of entropy 

For SL drift-diffusion models, systematic derivation procedures do exist. Many 
authors have used the Einstein relation to write a field-dependent diffusion co- 
efficient from the electron drift velocity. The latter was obtained by Ignatov and 
Shashkin 8 from a simple solution of a Bhatnagar-Gross-Krook (BGK) kinetic equa- 
tion 9 . Recently, we have used the Chapman-Enskog method to derive a general- 
ized drift-diffusion equation (GDDE) from a BGK-Poisson system 10 . Interestingly 
the Einstein relation does not hold, except in a very particular high-temperature 
limit 10 . Earlier, Ccrcignani, Gamba and Lcvcrmore had used essentially the same 
method to obtain reduced balance equations for a semiclassical BGK-Poisson semi- 
conductor system treating the band dispersion relation in the parabolic approxima- 
tion 11 . The aim of the present paper is to extend these methods to quantum kinetic 
equations. To this end, we shall derive a Wigner-Poisson system of equations for a 
SL having only one populated miniband, ignore changes in lateral momentum and 
model the collisions in the BGK approximation. The resulting quantum kinetic 
equation is 
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— e < are the one-particle Wigner function, the 2D electron density, the 2D doping 
density, the miniband dispersion relation, the barrier width, the well width, the SL 
period, the electric potential, the SL permittivity, the effective mass of the electron, 
the Boltzmann constant, the lattice temperature, the energy broadening of the 
equilibrium distribution due to collisions, the frequency of the inelastic collisions 
responsible for energy relaxation, the frequency of the elastic impurity collisions 
and the electron charge, respectively. The chemical potential \i is a function of n 
resulting from solving equation (1.3) with the integral of the collision-broadened 3D 
Fermi-Dirac distribution over the lateral components of the wave vector (k, k y , k z ), 
which is given by Eq. (1.4). Notice that, following Ignatov and Shashkin 8 , we have 
not included the effects of the electric potential in our Fermi-Dirac distribution. As 
discussed later, the model equations we use can be improved by including scattering 
processes with change of lateral momentum and an electric field-dependent local 
equilibrium. However the resulting model could only be treated numerically and 
the qualitative features of our derivation and nonlocal drift-diffusion equation would 
be lost in longer formulas. 

The Wigner-Poisson-Bhatnagar-Gross-Krook (WPBGK) system (1.1) to (1.4) 
should be solved for a Wigner function which is 27r/Z-periodic in k and satisfies 
appropriate initial and boundary conditions. If we integrate (1.1) over k and use 
the periodicity condition, only the first two terms in its left hand side yield nonzero 
contributions which, in fact provide the charge continuity equation as we shall see 
later. Equation (1.1) contains two pseudo-differential operators, depends on the 
unknown electric potential W and is nonlocal in x. However, the fact that / is 
2-7r/Z-periodic in k makes it possible to use a modification of the Chapman-Enskog 
method along the lines of that used in References 11 and 10 for the scmiclassical 
case and earlier in 12 and 13 for nonlinear Fokkcr-Planck equations. The result is a 
generalized drift-diffusion equation for the electric field which is nonlocal in space 
because it contains averages of the electric field over one or more SL periods. The 
nonlocality of the quantum drift- diffusion equation (QDDE) is a direct consequence 
of the nonlocality of the Wigner equation (1.1). 

In this paper, we derive the WPBGK system, deduce the nonlocal QDDE by 
means of the Chapman-Enskog method and solve it numerically to illustrate the 
resulting self-sustained oscillations of the current through a finite superlattice sub- 
ject to an appropriate dc voltage bias. The rest of the paper is organized as follows. 
In Section 2 we derive the WPBGK system. In Section 3, we study the hyperbolic 
limit of the WPBGK system and use the Chapman-Enskog method to derive the 
reduced equation to leading order. This equation is hyperbolic, and we would like 
to regularize possible singularities. Diffusive terms are obtained in Section 3 by 
calculating the first order terms in the Chapman-Enskog expansion. The resulting 
QDDE is discussed in Section 4 and solved numerically for appropriate voltage bias 
conditions. Finally, Section 5 contains a discussion of our results. 

2. Derivation of the quantum kinetic equation 
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Let us consider an n-doped SL formed by periodic alternation of two different semi- 
conductors such as GaAs (quantum well W) and AlAs (barrier B). The SL cross 
section S is much larger than I 2 , the square of the SL period. For n-doped SLs, we 
can restrict ourselves to studying electronic transport in the conduction band of the 
SL. We shall assume that the SL is under a dc voltage bias, which is equivalent to 
an external electric field directed along the SL growth direction. The corresponding 
Hamiltonian is 

H = Ho + H e - e + H 8C . (2-1) 

We have separated the electron-electron interaction H e _ e and other scattering pro- 
cesses (impurity, phonon,...) H sc from the one-electron Hamiltonian Ho- Typically, 
the electron-electron interaction is treated in the Hartree approximation. Then we 
can find the spectrum of the Hamiltonian H + H e - e by solving a non-linear station- 
ary Schrodinger-Poisson system of equations. Their solutions yield a basis in which 
quantum kinetic equations describing the scattering processes out of equilibrium 
can be written, as shown below. The envelope wave function is 14 

p„(x,k) = <p v (x,x±,k,k±) = -^=e iki - x ^„(x,fc,k ± ). (2.2) 



At zero external field, <p v satisfies 
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Here £ c is the conduction band edge of material W (GaAs, well), W(x) is the 
electric potential due to the electron-electron interaction, N 3 d{x) = Nd(x)/1 is the 
3D doping density, and 

, . _ J mw if x corresponds to a quantum well W . , 

\ tub H x corresponds to a barrier B 

. . J ew if x corresponds to a quantum well W ,_ 
^ ' 1 sb if a? corresponds to a barrier B 

are the masses and permittivities of the well and barrier. If V c corresponds to the 
conduction-band offset between the well and barrier material, we have 

V (x) — { ^ ^ X corres P on( ^ s t° a quantum well W ^ j\ 

c ^ [ V c if x corresponds to a barrier B. 

Moreover, the 3D equilibrium electron density is 

n {x) = \<f»(x,k,k ± )\ 2 n F (v,k,k ± ), (2.8) 
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where np denotes the Fermi function of the miniband v. The boundary conditions 
at the well-barrier interfaces are that ip„ and m(x)~ 1 dif v /dx are both continuous. 
The electronic spectrum is continuous, consisting of minibands of energies £ v {k; kj_) 
with v = 1,2,3, ... (doubly degenerate because of spin) with an associated basis of 
spatially extended Bloch wave functions <p v (x, k; kj_) = e u v (x, k; kj_) 14 . 

Although we can discuss the effects of scattering using the previous basis of 
electronic states that solve a non-linear Schrodingcr-Poisson system, we will for 
simplicity ignore the difference in mass and permittivity between barriers and wells 
and assume that the doping density is uniform. Of course in real devices, only the 
central parts of the wells are doped and the effective masses are different in wells and 
barriers 15 . These differences can be treated with the formulation sketched here (see 
for example, the Galdrikian-Birnir work for GaAs quantum wells, AlGaAs barriers 
and a nonuniform doping density 16 ), but then we would have to rely more heavily 
in numerical solutions than we intend to do, while obtaning qualitatively similar 
results (the AlAs barriers are much larger than the AlGaAs barriers and therefore 
the change in the energy levels due to electrostatic effects is much less important 
in the GaAs/ AlAs SLs we consider in this paper than in Galdrikian and Birnir's 
work). With these simplifications, the model [equations (2.3)-(2.7)] becomes the 
Kronig-Penney model with mw — mB = m* > £w = £b = £, and N^d^x) — N^d 
(constant). The Schrodinger-Poisson problem is simply [Ho — £,y(k)] y„(x, k) = 
with n n — N 3 d, and its solutions have the form 

^(x,k) = -i=e ik — <p v {x,k), £„(k) = ^+£,(fc). (2.9) 
The Bloch functions are 2-7r//-periodic in satisfying the orthogonality condition 

/oo 
<^(x, k) <p„{x, k') dx = 6^5(k - k'), (2.10) 
-OO 

and the closure condition 

/oo 
(Pp(x,k)tp v (x\k)dk = 8^5{x-x r ), (2.11) 
-oo 

provided the integral of \ifv\ 2 over one SL period is unity. 

Scattering different from electron-electron scattering is usually treated by writ- 
ing equations for the density matrix, its Wigner transform, or the non-equilibrium 
Green's function (NGF) i 7 ' 18 ' 1 . Whatever the chosen formulation, the equations 
for the one-electron functions depend on two-electron and higher functions, and 
we have the usual infinite hierarchy of coupled equations, which is well-known in 
classical kinetic theory. Typically, the hierarchy is closed by assuming some depen- 
dence of the two-electron functions on one-electron functions, which is suggested by 
perturbation theory in the limit of weak scattering 18 . Assuming weak scattering, 
the differences between the equations corresponding to the different formulations 
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are small. The trouble is that the kinetic equations are often used in the opposite 
hydrodynamic limit, in which collisions due to scattering are dominant. Then the 
results of using different formalisms are not equivalent, which has resulted in some 
discussion and confusion. In this paper, we shall not discuss the difference between 
formulations in a precise way. Instead, we shall write kinetic equations for the one- 
electron density matrix or Wigner function leaving unspecified the collision terms 
as much as we can, and discuss how to obtain reduced theories for electric field, 
electron density and current, and so on. These theories are easier to analyze and to 
solve numerically, and they are the ones commonly used to understand non-linear 
phenomena in SLs. 

To find a kinetic equation, we start writing equations for the coefficients a V! k(t) 
in the expansion of the wave function 



<A(x, t) = a*,k(t) ¥V(x, k) = M*, *)■ ( 2 - 12 ) 

If we ignore the scattering term H sc in equation (2.1), the coefficients a v ^(t) become 
d 
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Notice that equation (2.2) implies 

$„( x ,x') = S(x± -x'±)(() v (x,x'), 
<j> v {x, x') = ]T ip v {x, k) tpitf, k), (2.17) 

k 

and the closure condition in equation (2.11) yields 

^$„(x,x') = <Kx-x'). (2.18) 

Thus <&„(x, x') can be considered as the projection of the delta function <5(x — x') 
onto the band v. 
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After second quantization, the band density matrix is denned by 

p^(x,y,t) - (^(x,i)Vv(y,i)>, (2.19) 
so that the 2D electron density is (the factor 2 is due to spin degeneracy) 

n(x, t) - 21 5^j,(x, t) ^(x, t)> = 2? ^ p„ iW (x, x, t). (2.20) 

Using equations (2.19) and (2.20), we can derive the following evolution equation 
for the band density matrix 

d 2 d 2 



d_ _r 



m— — oo 

+eJ2 J W(z) [*„(y,z)p MV /(x,z,t)-$ M (z,x)^ 1 ,(z ) y ) t)]=Q[p] ) (2.21) 

with = in the absence of scattering. The Hartree potential satisfies the 
Poisson equation 



d 2 W e. 
e -mr = T (n - N D ) 



dx2 (2.22) 

(recall that e > and that W is the electric potential. The charge density is thus 
equal to minus the right hand side of this equation) . When considering scattering, 
the right hand side of equation (2.21) is equal to a non-zero functional of the band 
density matrix Q[p], whose form depends on the closure assumption we have made 
to close the density matrix hierarchy In the semiclassical limit, the kernel of the 
collision term Q[p] is usually found by using perturbation theory in the impurity 
potential, electron-phonon interaction, etc. For the time being, we shall not try to 
formulate collision models. Instead and in order to make contact with the kinetic 
equations in the semiclassical limit, we shall rewrite equation (2.21) in terms of the 
band Wigner function 

to M)1/ (x,k,t) = J Pf,., 

The evolution equation for the Wigner function is 
d Tiki d 
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in which the collision term is again left unspecified. Notice that the 2D electron 
density is 

11 f 

n(x, t) = / k, t) dk, (2.25) 

because of equation (2.20) and the definition in equation (2.23). From equations 
(2.24) and (2.25), we obtain the charge continuity equation 



<9J ie ^ 
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E v (m)w^ v [ x + ^,x_i_,k,i 
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dk, (2.28) 



provided our collision model satisfies / ^ ^ Qn,vdk = 0. 

A related formulation of the band Wigner functions (without collision terms) 
is due to Dcmcio et al 19 . One difficulty with our formulation is that the Wigner 
function in equation (2.23) is not 27r//-periodic in k. This can be corrected by using 
the following definition 

oo 
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together with the definition of equation (2.23). From equations (2.25) and (2.29), 
we obtain the 2D electron density in terms of /^^ 

n(x,t) = n V / / /„ lV (x,k,t)dfcdk ± . (2.31) 



Similarly, the transversal current density can be obtained from equations (2.27) and 
(2.29) 



J? / ^ £ ; ^ (x ' k ' * } * dk± - (2 - 32) 
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The current density along the growth direction has the form 
2eh — /■'/' 



J = 



87r 3 r 



fT/i r / 2irs\ ( 2lTS \ 
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from which we can also derive equation (2.28). 

The definition of periodic band Wigncr function is related to that adopted by 
Bechouche et al 20 . These authors have rigorously proved that the collisionless 
Wigner-Poisson equations for a crystal become the crystal Vlasov-Poisson equations 
in the semiclassical limit assuming that the initial conditions are concentrated in 
isolated bands. Scattering other than electron-electron scattering is not considered 
in these works. A similar work for a stratified material can be found in 21 . 

To find a Wigner-Poisson description of transport in a single miniband, we sum 
all the Wigner equations (2.24) over the band indices and use the closure condition 
in equation (2.18), so as to find an equation for w(x, k, t) = ^ v w M ,„(x, k, t) 



d_ 

dt 



Tiki 



w - 



d i 
dx± h 



e imkl 



m— — oo 



ml 



+ 



ie 



E v {m)w^ v ( x + ^,x_L,k, t 



'-27^- x - 



-W x 



2idk 



X_L 



w 



(2.34) 



Let us now assume that only the first miniband is populated and that there are 
no transitions between minibands, io(x, k, t) w wi.i(x, k, t). This approximation is 
commonly used when describing strongly coupled SLs with wide minibands. Then 
equation (2.34) becomes 
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This yields the following equation for the periodic Wigner function in equation 
(2.29) 



d Kk ± 

dt m* dx± 



e lmkl E x (m) 



m— — OO 



fix 



ml 



,xj_,k,t 



-fix 



ml 



,x_L,k,i 



ie 






W ^ 


+ H 



2i dk 



f = Q[f}- 



(2.36) 



10 Wigner- Poisson and nonlocal drift- diffusion model equations for semiconductor superlattices 



The dispersion relation S\(k) is an even periodic function of k with period 2tt/1 that 
can be written as £\{k) = Ai[l — cos(fc/)]/2 plus a constant in the tight-binding 
approximation (Ai denotes the width of the first miniband). Moreover, the field 
F = dW/dx (note that the real electric field is — F) satisfies 

e (S + d ! r Fi ) = T ( "- AW - <2 - 37) 

fOO poo 



I rn/l roc rOO 

n(a;,x_i_,i) = — — / / / /(a;,xj_, k,k.±,t)dkdk±. (2.38) 

^ J-TT/lJ-ooJ-oo 



We want to explicitly derive reduced balance equations from the kinetic equa- 
tion. For this purpose, we need a sufficiently simplified description of scattering. 
Scattering processes such as phonon scattering change the energy and momentum 
of the electrons leading the distribution function toward thermal equilibrium. We 
can describe these processes by a BGK collision model 9 

Qen[f] = -Ven(f-f FD ), (2.39) 

Jo [E - £ 1 (k) - ] + 1 + cxp (f^J 
n(x ) t) = - T / / / f FD (k,k ± ;n)dkdk ± . (2.41) 

471" J-TT/lJ-ooJ-oo 



Here T measures the finite width of the spectral function in thermal equilibrium 
due to scattering 1 . As T — > 0, the first factor in equation (2.40) becomes a delta 
function, and we recover the usual Fcrmi-Dirac distribution function with a chemical 
potential /x. The chemical potential /i = /j,(x,x±,t) is a function of the exact 
electron density n of equation (2.38) that is calculated by solving equation (2.41). 
With these definitions, the integral of Q e n[f] over momentum vanishes, and the 
equation of charge continuity holds, as we shall show below. Notice that we have 
not included the electric potential in (2.40). Then the equilibrium Wigner function 
does not include the h corrections to the semiclassical Fermi-Dirac distribution 
corresponding to the effective Hamiltonian £i(—id/dx) — (?i 2 /2m*)9 2 /9x^ — eW(x) 
(recall that the electron charge is — e): these corrections vanish if we set W = 
6,22 . Omitting the electric potential W is certainly an imperfection of our BGK 
collision model, which could perhaps be corrected using field-dependent collision 
models as in 6 . However, the resulting technical complications would encumber 
our derivation of drift-diffusion type equations. Thus we prefer to adopt the same 
line as in previous works 8 ' 10 and leave the study of field-dependent and broadened 
collision models for the future. 

Other processes, such as impurity scattering, conserve the energy of the electron, 
change only its momentum, and also preserve charge continuity. Gerhardts 23 used 
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the following model 

~ pir/l poo poo 

Qim P [f] = J J J[E(kM-E(k',k' ± )] 

x [f(k,k ± )-f(k',k' ± )]dk'dk' ± , (2.42) 

which can be rewritten as a relaxation to a weighted, energy-conserving mean 
value of the Wigner function: Qi mp [/] = —vg (/ — ®f[E(k)]), provided $f[E] = 
JS[E-E(k')}f(k')d 3 k'/[4n 3 g(E)}, g(E) = J 5[E - E(k')} d 3 k'/(4n 3 ) is the density 
of states, and vq = Vi mp g(E(k)). Note that we have dropped the dependence of the 
Wigner function on space and time. This collision term couples the vertical mo- 
tion of the electron to the lateral degrees of freedom. For SLs, scattering processes 
involving acoustic phonons, impurities and interface roughness may modify the lat- 
eral momentum. Liu et al 24 have calculated and measured the scattering times 
in weakly coupled double barrier heterostructures and obtained energies smaller 
than 10 meV. This energy is also a generous upper bound for the energy exchanged 
in impurity scattering processes in strongly coupled SLs, as indicated in Wacker's 
review (his smallest scattering time is 0.0666 ps 1 ). In any case, 10 meV is much 
less than the typical energy exchange in the direction parallel to the electric field, 
of more than 1 eV (which can be estimated from the width of the Brillouin zone, 
2w/l). Thus, we may reasonably assume that the variation of the energy in the 
lateral direction is negligible, E(k,k±) — E(k',k' ± ) rts £\{k) — £i(k'), and therefore 

Qimp[/] « SMk) [f(k) f(k')} dk' 

= [m s ^ [m fi ~ k)] ' (2 ' 43) 

f( x ' k > t ) = ^2s J J f{x,x±,ky ± ,t)dk! ± dx ± . (2.44) 

Here £[{k) — d£\jdk. The approximate collision term in equation (2.43) was intro- 
duced by Ktitorov et al 25 . In terms of the ID Wigner function in equation (2.44), we 
can derive from equations (2.36)-(2.44) the ID WPBGK system (1.1)-(1.4) (which 
is independent of xj_ if we assume that the initial Wigner function does not depend 
on xj_). The Gerhardts collision model is still better than the approximation (2.43) 
because his collision model explains the observed temperature dependence of the 
field at which the drift velocity has its lowest local maximum 23 . However the sim- 
pler ID WPBGK model is much more convenient for analytical calculations. The 
scmiclassical limit of the WPBGK system has been analyzed in 10 in order to derive 
a reduced drift-diffusion model in the hydrodynamic limit. Particular solutions of 
the semiclassical model have been found by different authors. Ignatov and Shashkin 
8 found the stationary, space- independent semiclassical solution of equation (1.1) 
for a Boltzmann distribution function, the Boltzmann limit of equation (1.4), and 
studied its linear stability to plane wave disturbances in the field. 
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Before we proceed, it is convenient to derive the charge continuity equation and a 
form of Ampere's law for the current density from equation (1.1). Since the Wigner 
function is periodic in k, we can write the second and third terms on the right hand 
side of equation (1.1) in terms of its Fourier series 



f(x,k,t)= fj(x,t)e iikl - 



(2.45) 



We have 



,ijkl 



(2.46) 



and therefore, 



W x 



_L__9 

2idk' 



t)f-W(x 



)_d_ 

2idk : 



,t)f = 



kl 



kl 



Y, [W(x + jl/2, t) - W(x - jl/2, t)} / je « 

j=-oo 

oo 
j= — oo 

Here we have defined the average 

(F)j(x, t) = - / F{x + s,t)ds. 

.]<■ J -jl/2 

Note that differentiating an average, we obtain a finite difference 

d_, \ _.dg, _ g(x + jl/2) -g(x- jl/2) 
dx {9)j ~ [ dx )j ~ 



Jl 



Then the second term in equation (1.1) is 



Jl 



f [x + — ,k,t) - f [x- —,k,t 



- y 

j = -oo 



Jl 



j = -oo 



(2.47) 



(2.48) 



(2.49) 



(2.50) 



which in the case of the tight-binding dispersion relation £\(k) — Ai (1 — coskl)/2 
becomes v(k) d(f)i/dx, with the usual miniband group velocity 



(2.51) 
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Inserting equations (2.47) and (2.50) into equation (1.1), we obtain the following 
equivalent form of the Wigner equation, which is particularly suitable for treating 
SL problems 



§£+£; f e «( £l0) | (/)j+e(f)j/j ) 



j = -oo 

= -u m (/ - f D ) - * mp /^iM^iZM , (2 .52) 

We now integrate this equation over k, thereby getting the charge continuity equa- 
tion 

We can eliminate the electron density from equation (2.53) by using the Poisson 
equation (1.2) and integrating the result over x, thereby obtaining the non-local 
Ampere's law 

£ % + 1 jpxMEii-mh = m- (2-54) 

Here J(t) is the total current density. Equations (2.52), (2.53) and (2.54) are spa- 
tially non-local versions of the corresponding semiclassical equations. The charge 
continuity and Ampere's equations have their traditional form as derived from semi- 
classical Boltzmann equations, except that the electron current is averaged over the 
SL periods. This non-locality will be transmitted to the QDDE. 

3. Derivation of the QDDE in the hyperbolic limit of the kinetic equation 

To derive the QDDE, we shall assume that the electric field contribution in Eq. 
(2.52) is comparable to the collision terms and that they dominate the other terms 
(the hyperbolic limit) 10 . Let vm and Fm be electron velocity and field scales 
typical of the macroscopic phenomena described by the sought balance equation; 
for example, let them be the positive values at which the (zeroth order) drift 
velocity reaches its maximum. In the hyperbolic limit, the time to it takes an 
electron with speed vm to traverse a distance x — eFmI/(bNd), over which 
the field variation is of order Fm, is much longer than the mean free time be- 
tween collisions, ~ %/(eFMl) — t\. We therefore define the small parameter 
A = h/to — HvmNd/(sF m 1 2 ) and formally multiply the first two terms on the left 
side of (1.1) and of (2.52) by A 10 . After obtaining the number of desired terms, we 
set A = 1. 

The solution of equation (2.52) for A = is the stationary space-independent 
solution that is easily found as a Fourier series 

/(°)(fc;F)= £ /fe«« /f = \"^g4 Te /r, (3-1) 

j=-oo 3 
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in which 



F M 



Tt^/Ven{Ven + ^imp) 



Ven + Vu 



(3.2) 



Since f FD is an even function of k, its Fourier coefficient f? D is real. Note that 
equation (1.3) implies = f FD = n. 

The Chapman-Enskog ansatz consists of writing the distribution function as an 
expansion in powers of the book-keeping parameter A (recall that we have to set 
A = 1 after retaining the desired number of terms) 10 



/(*, k, t; A) = /(°)(fc; F) + J2 f (m) F ) 

of 



m— 1 



dt 



J2 J (m) {F)\ m = Jit). 



(3.3) 
(3.4) 



m=0 



The coefficients f^ m ^{k;F) depend on the 'slow variables' x and t only through 
their dependence on the electric field and the electron density (which are related 
through the Poisson equation). The field obeys a reduced evolution equation (3.4), 
in which the functionals j( m ) (F) are chosen so that the /( m ) (k; F) are bounded and 
27r/Z-periodic in k. Differentiating Ampere's law (3.4) with respect to x, we obtain 
the charge continuity equation. Moreover the condition 



/ f( m \k;n)dk = 2irf( m) /l = 0, m > 1, 

J-tt/I 



(3.5) 



ensures that /( m ) for m > 1 does not contain contributions proportional to the 
zero-order term f^. Note that the insertion of equation (3.3) in Ampere's law 
(2.54) yields 



JM = £ £ j<Ml*(-j)/, (m) ]>j, 



h 



(3.6) 



which is also obtained by means of the above mentioned boundedncss condition. 

Inserting equations (3.3) and (3.4) into equation (2.52) and equating all coeffi- 
cients of A m in the resulting series to zero, we find the hierarchy 



£/ (D 



(2) : 



dt 
dt 



— -co 



d 



(0)\ 



lf(0) 

at'' 



(3.7) 
(3.8) 
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and so on. We have defined 

Lu{k) s * ± jl(F) jUj e« kl + (v en + ^p) u(k) - ^- k \ (3.9) 

j=-oc 

and the subscripts and 1 in the right hand side of these equations mean that 
edF/dt is replaced by J - j(°\F) and by —J ( - 1 \F), respectively. 

The linear equation Cu = S has a bounded 27r/l-periodic solution provided 
J^J^Sdk = 0. Equation (3.7) and this solvability condition yield equation (3.6) 
for m = 0. The solution of equation (3.7) is 

m ~ ReS^ +iT- 2 lmS ( j 1) - ijfjS^ /r e .... 
fV(k;F)= ]T —1 > 1 > e e«", (3.10) 

in which S*^ is the jth Fourier coefficient of the right hand side of equation (3.7). 
Using equation (3.10), we can now explicitly write two terms in equation (3.4), 
thereby obtaining the following QDDE for the field and the electron density given 
by the Poisson equation (1.2) (with d 2 W/dx 2 = dF/dx). For the tight-binding 
dispersion relation, the QDDE is 

N = (nVM)i + ((A - l){{nVM)i)i)i 
ZT e Ai / B 



F M K{v en + v imv ) \l + ^ 2 /i' (3 ' 13) 

D _ l 2 ^l ( d 2 (F)i ^v M r e C \ 

8h 2 (v en + is imp )(l + T 2 ) \ dx 2 I A, ' 



n , ±T 2 nM 2 d(F) 2 
" ~~ (1+4^|) 2 8x / 1 

nM 2 (l-4T 2 )d(F) 2 



+I s (1 + 4^|) 2 dx , 
4hv M (l + T 2 )F(nMy I .. \-T 2 0(F) 



/T e (l + J- 2 )A! \ {l + T 2 ) 2 dx 

/ (nM 2 )' 8 2 F\ / (nM 2 )^2 d 2 F 

\ 1 + Ajq dx 2 I x \ 1 + 4JT| dx 2 



nM - , -rr. \' L (3.16) 
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8hv M {l + T^){nMy T I (nM)'T d 2 F \ 
+ /r e (l+JT2)A 1 \ l + T 2 dx 2 / 1 



(3.17) 



mmn D ) = ™*m Mi{n/ND)= *mm t (3 . 18) 




la (E-S + Scosk) 2 + T 2 



dk. (3.19) 



Here g' denotes dg/dn, S = A 1 /(2k B T), p, = ^/{k B T), f = T/(k B T), T = T x , and 
n = Nd at the particular value of the dimensionlcss chemical potential ft = M. If 
the electric field and the electron density do not change appreciably over two SL 
periods, (F)j w F, the spatial averages can be ignored, and the non-local QDDE 
(3.11) becomes the local generalized DDE (GDDE) obtained from the semiclassical 
theory 10 . The boundary conditions for the QDDE (3.11), which contains triple 
spatial averages, need to be specified for the intervals [—21, 0] and [Nl, Nl + 21], and 
not just at the points x = and x = Nl (N denotes the number of SL periods 
spanning the device), as in the case of the parabolic semiclassical GDDE. Similarly, 
the initial condition has to be defined on the extended interval [—21, Nl + 21]. Note 
that the spatial averages in the nonlocal QDDE give rise to finite differences of 
partial derivatives in the diffusion terms, and therefore lead to a type of equations 
for which little seems to be known. 



4. Numerical solution of the QDDE 

In this Section, we solve the nonlocal QDDE (3.11) together with the voltage bias 
condition 

pNl 

I F(x,t)dx = 4>F M Nl, (4.1) 
Jo 

for the field F and the total current density J. As boundary conditions in the 
intervals [-21,0] and [Nl,Nl + 21], we adopt 

dF 

Js^=aF, (4.2) 

at all points [ — 2i, 0] of the Ohmic injecting contact and zero-flux boundary condi- 
tions at the receiving contact [Nl, Nl + 21]. The contact conductivity a is selected 
so that eN r>v mV (F / Fm) and aF intersect on the second branch of V(T), in which 
dVjdT < 0. This is a typical boundary condition yielding self-sustained oscilla- 
tions in drift-diffusion SL models 2 ' 10 . We have used a constant initial condition 
F = 4>Fm in our numerical simulations. The SL parameters we have used corre- 
spond to a 157-pcriod 3.64 nm GaAs/0.93 nm AlAs SL 26 at 14K, with A = 72 
meV, N D = 4.57 x 10 10 cm~ 2 , v imp — 2v en — 18 x 10 12 Hz under a dc voltage 
bias of 1.62 V (cfi = 1). Cathode and anode contact conductivities are 2.5 and 0.62 
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Figure 1: (a) Current (Jo = cvmNd/1) vs. time (to = eFm/Jq) during self- 
oscillations for a voltage biased GaAs/AlAs SL, as described by the QDDE (solid 
line), the QDDE with T = (long-dashed line) and by the GDDE (dot-dashed line), 
(b) Comparison between the fully developed dipole wave for the QDDE (solid line) 
and the dipole wave for the GDDE (dashed line) . (c) Dipole wave at different times 
during the stage in which it is shed from the injecting contact, (d) Same as (c) 
for the stage in which the dipole disappears at the anode, located at NI/xq « 44. 
Parameter values are xq = vm^o = 16 nm, to = 0.43 ps, Jo = 6.07 x 10 5 A/cm 2 , 

4>=i. 



Our numerical solution shows that the current and the field profile become 
stationary for <f> < 0.75 (1.2 V). For larger values of the dimensionless voltage (f>, 
the initial field profile evolves toward a stable time-periodic solution for which J 
oscillates with time and the field profile shows recycling and motion of a pulse from 
x = to the SL end. Fig. 1 shows the self-oscillations of the current for 1.62 
V ((f> = 1) and the corresponding field pulse at different times. In this figure, we 
compare the solution of the GDDE corresponding to the semiclassical BGK-Poisson 
kinetic equation and the solution of the QDDE for T = (no collision broadening) 
and for T = 18 meV which is of the same order as the collision frequencies. 

Self-oscillations in the QDDE have a frequency of vq = 25.5 Ghz, faster than in 
the GDDE, v G = 20.6 Ghz (relative frequency (vq-v g )/v g = 23.8%). When T = 0, 
the frequency is vq* = 21 Ghz, and the relative frequency (vq» — vg)/ v g = 1-94%. 

Collision broadening shortens the period of the current oscillations and therefore 
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it reinforces the effects of the nonlocal terms in the QDDE due to quantum effects. 
5. Discussion 

We have derived a nonlocal drift-diffusion model for a strongly coupled superlat- 
tice with one miniband by applying the Chapman-Enskog procedure to a simple 
WPBGK quantum kinetic equation. Our local equilibrium function is a collision- 
broadened Fermi-Dirac distribution inspired in the well-known thermal equilibrium 
of the Kadanoff-Baym equations 17,1 , but the effects of the electric field on the 
collision operator have been ignored. For the case of a quantum particle in an 
external potential, Degond and Ringhofcr 6 have proposed a different form of the 
BGK collision model, in which the density matrix corresponding to local equi- 
librium is a Maxwellian operator obtained by minimizing the Boltzmann entropy 
subject to moment constraints. Their formulation includes the effects of the elec- 
tric field on the collisions, and the chemical potential is a nonlocal functional of 
the electron density, not a function given by solving Eq. (1.3) for /x. This func- 
tional relation introduces additional nonlocalities in the balance equations resulting 
from the Chapman-Enskog procedure. See 27 , in which quantum drift-diffusion and 
energy-transport models are derived for the case of a particle in an external po- 
tential undergoing collisions according to their BGK model. Notice that the BGK 
collision models used by these authors conserve kinetic energy, which is a somewhat 
unrealistic feature for semiconductors. 

In this paper, we have treated the case of a one-dimensional SL with one mini- 
band. We have used a Wigner-Poisson equation with relaxation towards a broad- 
ened local Fermi-Dirac distribution and a simple elastic impurity collision term. We 
have derived a nonlocal drift-diffusion equation and solved it numerically to show 
that our model displays the self-oscillations of the current which are observed in 
experiments. There are several directions in which we could extend our results. 
Firstly, it would be interesting to derive quantum hydrodynamic models by using 
a more general BGK collision model and compare it with the models by Lei and 
collaborators 3 ' 4 ' 5 . Our approximation /(°) is possibly closer to the real electron dis- 
tribution than the parametrized local equilibrium used by Lei et al, which has been 
criticized by other authors 1 . Unlike Lei's, our approximation /(°) takes explicitly 
into account the effect of the electric field (and therefore the modified "equilibrium" 
f(0) correc ts somewhat the effect of having neglected the electric potential in our 
collision model) . Secondly, it would be important to consider the case of a SL with 
several minibands, although figuring out a reasonable local equilibrium distribution 
to include in the BGK collision model is probably harder. Thirdly, there are many 
interesting mathematical questions arising from these derivations. To mention one, 
the trailing edge of the pulse in Fig. 1 is rather sharp. In related work with drift- 
diffusion models of the Gunn effect, it is known that this trailing edge is close to a 
shock wave of the zero-order hyperbolic drift equation 28 . Obviously, the diffusive 
terms in the QDDE regularize this shock, and so do the collision terms in the kinetic 
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equation (even if we think in terms of the semiclassical kinetic equation) . What is 
the relation between these two regularizations? Cafiisch has considered a related 
problem in gas dynamics 29 . We hope that our work paves the way to tackling these 
problems in the future. 

Acknowledgment 

This work has been supported by the MCyT grant BFM2002-04127-C02-01, and by 
the European Union under grant HPRN-CT-2002-00282. R.E. has been supported 
by a postdoctoral grant awarded by the Consejerfa de Education of the Autonomous 
Region of Madrid. 

References 



1. A. Wacker, Semiconductor superlattices: A model system for nonlinear transport, Phys. 
Rep. 357 (2002), 1-111. 

2. L. L. Bonilla, Theory of Nonlinear Charge Transport, Wave Propagation and Self- 
oscillations in Semiconductor Superlattices, J. Phys.: Condens. Matter 14 (2002), R341- 
R381. 

3. X. L. Lei and C. S. Ting, Theory of nonlinear electron transport for solids in a strong 
electric field, Phys. Rev. B 30 (1984), 4809-4812. 

4. X. L. Lei, Distribution function and balance equations of drifting Bloch electrons in an 
electric field, Phys. Rev. B 51 (1995), 5526-5530. 

5. J. C. Cao and X.L. Lei, Synchronization and chaos in miniband semiconductor super- 
lattices, Phys. Rev. B 60 (1999), 1871-1878. 

6. P. Degond and C. Ringhofer, Quantum moment hydrodynamics and the entropy prin- 
ciple, J. Stat. Phys. 112 (2003), 587-628. 

7. K. Huang and X. Wu, Distribution function of a drifting electron gas for general energy- 
band structures, Phys. Rev. B 49 (1994), 2223-2225. 

8. A. A. Ignatov and V.I. Shashkin, Bloch oscillations of electrons and instability of space- 
charge waves in semiconductor superlattices, Sov. Phys. JETP 66 (1987), 526-530. 

9. P. L. Bhatnagar, E. P. Gross, and M. Krook, A Model for Collision Processes in Gases. 
I. Small Amplitude Processes in Charged and Neutral One-Component Systems, Phys. 
Rev. 94 (1954), 511-525. 

10. L.L. Bonilla, R. Escobedo and A. Perales, Generalized drift- diffusion model for miniband 
superlattices, Phys. Rev. B 68 (2003), 241304(R) (4 pages). 

11. C. Cercignani, I. M. Gamba and C. D. Levermore, A drift- collision balance for a 
Boltzmann-Poisson system in bounded domains, SI AM J. Appl. Math. 61 (2001), 1932- 
1958. 

12. L. L. Bonilla, Chapman- Enskog method and synchronization of globally coupled oscilla- 
tors, Phys. Rev. E 62 (2000), 4862-4868. 

13. L. L. Bonilla y J. Soler, High-field limit of the Vlasov-Poisson-Fokker-Planck system: 
A comparison of different perturbation methods, Math. Mod. Meth. Appl. Sci. (M 3 AS) 
11(8) (2001), 1457-1468. 

14. G. Bastard, Wave mechanics applied to semiconductor heterostructures. (Hal- 
sted, New York, 1988). 

15. H. T. Grahn, ed., Semiconductor Superlattices: Growth and Electronic Prop- 
erties, (World Scientific, Singapore, 1995). 



20 Wigner- Poisson and nonlocal drift- diffusion model equations for semiconductor superlattices 



16. B. Galdrikian and B. Birnir, Period doubling and strange attractors in quantum wells, 
Phys. Rev. Lett. 76 (1996), 3308-3311. 

17. L. P. Kadanoff and G. Baym, Quantum Statistical Mechanics: Green's func- 
tion methods in equilibrium and nonequilibrium problems. (Addison- Wesley, 
Redwood City, 1989). 

18. H. Haug and A.-P. Jauho, Quantum kinetics in transport and optics of semi- 
conductors. (Springer, Berlin, 1996). 

19. L. Demeio, L. Barletti, A. Bertoni, P. Bordone, C. Jacoboni, Wigner- function approach 
to multiband transport in semiconductors, Physica B 314 (2002), 104-107. 

20. P. Bechouche, N. Mauser and F. Poupaud, Semiclassical limit for the Schrddinger- 
Poisson equation in a crystal, Comm. Pure Appl. Math. 54 (2001), 851-890. 

21. P. Bechouche and F. Poupaud, Semiclassical limit of a Schrddinger equation for a 
stratified material, Monatsh. Math. 129 (2000), 281-301. 

22. C. L. Gardner, The quantum hydrodynamic model for semiconductor devices, SIAM J. 
Appl. Math. 54 (1994), 409-427. 

23. R. R. Gerhardts, Effect of elastic scattering on miniband transport in semiconductor 
superlattices, Phys. Rev. B 48 (1993), R9178-R9181. 

24. H. C. Liu, M. Buchanan, G.C. Aers, Z. R. Wasilewski, W.T. Moore, R.S.L. Devine and 
D. Landheer, Transport measurements of resonant tunneling widths, Phys. Rev. B 43 
(1991), 7086-7090. 

25. S.A. Ktitorov, G.S. Simin and V.Ya. Sindalovskii, Bragg reflections and the high- 
frequency condictivity of an electronic solid-state plasma, Sov. Phys. Solid State 13 
(1972), 1872-1874. 

26. E. Schomburg, T. Blomeier, K. Hofbeck, J. Grenzer, S. Brandl, I. Lingott, A. A. Ignatov, 
K. F. Renk, D. G. Pavelev, Y. Koschurinov, B. Y. Melzer, V.M. Ustinov, S. V. Ivanov, A. 
Zhukov, and P. S. Kopev, Current oscillations in superlattices with different miniband 
widths, Phys. Rev. B 58 (1998), 4035-4038. 

27. P. Degond, F. Mehats and C. Ringholer, Quantum hydrodynamic models derived from 
the entropy principle, Contemp. Math. , (2004), to appear. 

28. L. L. Bonilla, Solitary waves in semiconductors with finite geometry and the Gunn 
effect, SIAM Journal on Applied Mathematics 51 (1991), 727-747. 

29. R. E. Caflisch, Fluid dynamics and the Boltzmann equation, in Nonequilibrium phe- 
nomena I The Boltzmann equation ( J. L. Lebowitz and E. W. Montroll, eds) 
194-223. North Holland, Amsterdam 1983. 



